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ABSTRACT 

We present mid-infrared spectra and photometry of thirteen redshift 0.4 < z < 1 dust-reddened 
quasars obtained with Spitzer IRS and MIPS. We compare properties derived from their infrared spec- 
tral energy distributions (intrinsic AGN luminosity and far-infrared luminosity from star formation) 
to the host luminosities and morphologies from HST imaging, and black hole masses estimated from 
optical and/or near- infrared spectroscopy. Our results are broadly consistent with models in which 
most dust reddened quasars are an intermediate phase between a merger-driven starburst triggering a 
completely obscured AGN, and a normal, unreddened quasar. We find that many of our objects have 
high accretion rates, close to the Eddington limit. These objects tend to fall below the black hole mass 
- bulge luminosity relation as defined by local galaxies, whereas most of our low accretion rate objects 
are slightly above the local relation, as typical for normal quasars at these redshifts. Our observations 
are therefore most readily interpreted in a scenario in which galaxy stellar mass growth occurs first by 
about a factor of three in each merger /starburst event, followed sometime later by black hole growth 
by a similar amount. We do not, however, see any direct evidence for quasar feedback affecting star 
formation in our objects, for example in the form of a relationship between accretion rate and star 
formation. Five of our objects, however, do show evidence for outflows in the [O III]5007A emission 
line profile, suggesting that the quasar activity is driving thermal winds in at least some members of 
our sample. 

Subject headings: Quasars: general — infrared: galaxies — galaxies: starburst — galaxies: evolution 



1. INTRODUCTION 

The study of quasars across cosmic time has proven 
to be important as their evolution and growth is in- 
timately linked to that of their host galaxy. Mergers 
have long been invoked as the primary ignition mech- 
anism of q uasar activity, at least for the high luminos- 
ity objects dSanders fc MirabellU996l: lUrrutia , et"ai]|2008t 
iBennert et all 120081) . The loss of angular momentum 
during a gas rich merger allows the gas to be funneled 
to the center igniting a starburst and fueling the central 
black hole. In this model, eventually the quasar grows 
strong enough and develops a feedback mechanism that 
expels the obscuring material and shuts off star formation 
(jSilk fc Reeslll998tlHopkins et alll2008[ ). This scenario is 
further supported by the existence of tight correlations 
between the properties of galaxy bu l ges and their central 
black holes (Ma gorrian et all 119981: iFerrarese fc Merrittl 
120001) . 

Putting the picture together observationally is chal- 
lenging, however. Populations of quasars selected via dif- 
ferent techniques, although overlapping, differ substan- 
tially in their detailed properties and may bias towards a 
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particular type of quasar or evolutionary phase e.g. AGN 
selected in the infrared may belong to different class or 
evolutionary p hase than AGN sel ected in the radio or 
optical regime (jHickox et all 120091) . In an ideal survey 
we would be able to select quasars at specific points in 
their lifetimes. For example, to test a starburst- AGN 
connection we would prefer to select young, recently ig- 
nited objects to see if there is still star formation occur- 
ring in the host galaxies, and to identify likely feedback 
mechanisms. 

One approach to finding candidates for young quasars 
is to search for dust obscured quasars that may cor- 
respond to objects yet to fully clear out the dust and 
gas surrounding them. Conventionally, obscured quasars 
are divided into two classes: (1) the so-called type- 2 
quasars - objects showing only narrow emission lines 
with large inferred extinctions (Ay > 5 — 100) towards 
the nucleus, most likely caused by a torus surround- 
ing the accretion disk as require d in the AGN "unifica- 
tion by orientation" mo del (e.g.. lAntonucci et al.lfT993l : 
lUrrv fc PadovanH [1995) and (2) moderately reddened 
quasars (Ay ~ 1 — 5), which still show broad emis- 
sion lines in the rest-frame optical, and whose contin- 
uum is still dominated by the quasar rather than the 
host galaxy. These latter objects are obscured by a cold 
absorber along the line of sight to the quasar, most likely 
in the host galaxy. They may thus represent the young 
objects in the final stages of emerging from their dusty 
cocoons. Throughout this paper the terms red, dust- 
obscured and moderately reddened quasars all refer to 
the latter category of objects. 

Dust obscured quasars of all types are missing from, or 
are severely underrepresen ted in optical surv eys, but are 
present in the radio (e.g., iWhite et al.l 12003). in the X- 
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ray (e.g.. iHickox et aL"l l2007l). when sele cted on the basis 
of na rrow optical AGN emission lines (TZakamsk a et al 




, or when selected in the near-infrared ( Cutri et al 



Mid-infrar e d sur v eys (e.g.. lLacv et all 



iStern et al.1 120051 : lAlonso-Herrero et al 



2004, 



2006 



Donlev et al.l 120071 ) have also been extremely successful 



at finding and identifying obscured quasars over a wide 
range of reddenings, redshifts, and luminosities. Sur- 
veys combining deep radio and mid-infrared data show 
that the optically- and/or X-ray-selected quasar popu- 
lation const itutes less than half of the total population 
of qu asars (|Donlev et al.l 120051 : iMartmez-Sansigre et al.l 
2005). Joint selection using radio and near-infrared, as 
used for the sample in this paper, has proven to be one of 
the most reliable ways to select the moderately obscured 
quasar population. The requirement of a bright radio 
source and very red optical through near-infrared colors 
results in a set of candidates with relatively little contam- 
inatio n from normal ga l axies, low luminosity AGN and 
stars (iGregg et al Il2002t lLacv et al.ll2002t iGlikman et al.l 
l200l 120071 120121: lUrrutia et al.ll2009fl .~ 

The success of the Spitzer Space Telescope 
(| Werner et al.1 120031 ) has allowed for the detailed 
study of quasars in the mid-inf rared. Programs w ith 
the Infrared Spectrograph (IRS; iHouck et all 120041) . in 
particular, have refined our knowledge of QSO spectral 
energy distributions (SEDs) in the mid-infrared. For 
example, the detection of silicate in emission in IRS 
spectra of quasars i s a strong support for the un if icatio n 
model (jHao et al.1 120051: ISiebenmorgen et al.1 [2005). 
Average prop erties of classes of objects using IRS 
spectra (e.g, iSpoon et al.l I2007D have allowed us to 
constrain the physi cal properties of dust in AGN (e.g., 
iNikutta et all[200l . On average the spectra of luminous 
quasars are flat, show little or no PAH emission and the 
silicate features are in emission, associated with dust 
re-emission. In contrast, the spectra of Ultraluminous 
Infrared Galaxies (ULIRGs) show a steep rise towards 
the long wavelengths, moderate PAH emission and 
silicate absorpt ion troughs asso ciated with embedded 
star formation (jHao et al.l 12007). However, some type 
2 and some reddened quasars at moderate t o high red- 
shifts s how deep silicate absorpt i on fea t ures (lLacv et al 



2007bi IMartmez-Sansigre et al.l 120081: iZakamska eta!. 
20081) in contrast to samples of X-ray sele cted Type- 2 



quasa rs, in which they are generally weak ([Sturm et all 
I20M) . 

The Multiband Imaging Photometer for SIRTF (MIPS; 
iRieke et "al"1l2004D on board Spitzer has also allowed us to 
sample the colder part of the galaxy SEP . While AGN 
domi nate the bright 24/im population ([Donlev et al.l 
2008), the Spitzer 70/zm field population is domi- 
nated by ULIRGs whic h have a high merger fraction 
([Kartaltepe et al.ll2010al[ bT). However, AGN and quasars 
still make up a significant fraction of 70/im sources, and 
high star formation rates in quasar hosts are common 
(e.g. iSchweitzer et al.ll20"06l: iNetzer et al.ll2007l) . Consis- 
tent with those results iShi et al.l (|2009t )~find that Type-1 
quasars show star formation rates higher than the field 
galaxies, with luminosities typical of Luminous Infrared 
Galaxies (LIRGs, Lm — 10 11_12 L Q ). In particular, dust 
reddened quasars have been found to have higher than 
usual 60/12/^m lumi nosity ratios, an indicato r of higher 
star formation rates (Gcorgak akis et aLll2009D . 



In this paper we will describe Spitzer IRS and MIPS 
observations of 13 redshift 0.4 < z < 1 reddened quasars 
selected using joint radio and near-infrared selection for 
which we have also obtained Hubble S pace Telescope 
(HST) observations ([Urrutia et al.l 120081 ). We use these 
observations to estimate the intrinsic luminosities of the 
quasars, and to estimate or constrain the star formation 
rates in the host galaxies. We then evaluate the evidence 
that this population corresponds to an intermediate stage 
in quasar evolution by estimating black hole masses and 
accretion rates, comparing our HST host galaxy lumi- 
nosities and morphologies to the other properties, and 
by searching for evidence of feedback processes at work 
in the hosts. Throughout this paper we adopt a flat Uni- 
verse, Hq = 70 km s -1 Mpc -1 , Ha = 0.7 cosmology. 

2. QUASAR SAMPLE AND OBSERVATIONS 

In the past few years we have been selecting a sample 
of luminous lightly-obscured qu asars using a com bina- 
tion of the FIRST radio surv ey dBecker et al. l [i995l) . the 
2MASS near infrared survey ([Skrutskie et al.|[2006D and 
the Cambridge APM scans of the POSS plates. Using 
a J — K > 1.7 and R — K > 4 color wedge for the 
matched point sources, we have been very effective in 
finding these so-called red quasars with sp ectroscopic fol- 
low up at Keck (ESI) and IRTF (Spex) ([Glikman et al.l 
12004 120071: lUrrutia et all 120091) . Even though the sur- 
vey is radio selected, most of the objects rather fall in 
the radio intermediate regime than radio-loud, so it is 
unlikely that the redness of the objects comes from a 
synchrotron component, but is in fact associated with 
dust-reddening. Furthermore, the excess of steep spec- 
trum sources at faint radio fluxes and the point-source 
nature of the radio sources leads us to expect that the na- 
ture of the radio emission is due to a q uasar and not from 
star formation ([Glikman et al.l 120071 ) . As of early 2012, 
we have well over 120 spectrosc opically confirmed r ed 
quasars found with this method ([Glikman et al.ll2012f) . 

These red quasars are highly luminous and only moder- 
ately reddened (Ay ^1 — 5 mags) , so they outshine their 
host galaxies in the rest-frame optical, and are thus not 
red due to galaxy starlight. Their spectral energy dis- 
tributions (SEDs) are generally well fit by reddening a 
standard quasar template spectrum using an SMC dust- 
reddening law. It is only in the rest-frame UV, where 
the quasar light is heavily extinguished, that the host 
galaxy s tarts to play a role in t erms of overal l flux d en- 
sity. See IGlikman et al.l ([20071 ): lUrrutia et all ([2009D for 
in depth discussions of how reddening is estimated for 
these objects. 

A representative sample of 13 dust-reddened quasars 
with 0.4 < z < 1 were followed up with the ACS Wide 
Field Camera on HST with the F475W and F814W fil- 
ters. These are very luminous quasars at z ~ 0.7 hav- 
ing intrinsic absolute magnitudes in the range —23.5 > 
Mb > —26.2. The host galaxies show a high amount of 
interaction: 85% of the images show evide nce of mor- 
phological disturbance ([Urrutia et al.l [2008). We also 
observed that the more reddened the objects were, the 
more disturbed their morphologies (as measured by the 
Gini coefficient combined with the Con centration Index 
(|Lotz et all 120041: 1 Abraham et al.|[2003D V These results 
support the merger-induced origin for high luminosity 
quasars. 
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Table 1 

Spitzer IRS and MIPS observations 



Source IRS AORID MIPS AORID 524 (mjy) 5Vo (mJy) Sieo (mjy) 



F2M0729+3336 


22386176 


22386432 


12.9±3.0 


<15 


<30 


F2M0825+4716 


22389760 


22386688 


36.1±4.7 


85±19 


46±11 


F2M0830+3759 


22390016 


22386944 


26.1±3.9 


32±12 


<30 


F2M0834+3506 


22390272 


22387200 


17.6±3.4 


26±13 


<30 


F2M0841+3604 


22390528 


22390528 


11.7±3.0 


53±14 


<30 


F2M0915-2418 


17540352+ 


22387712 


87.3±7.0 


207±24 


73±12 


F2M1012+2825 


22390784 


22387968 


4.9±0.8 


<15 


<30 


F2M1113+1244 


22391040 


22388224 


59.4±5.8 


138±25 


67±12 


F2M1118-0033 


22391296 


22388480 


24.4±4.0 


80±24 


44±10 


F2M1151+5359 


22391552 


22388736 


6.9±0.7 


20±7 


<30 


F2M1507+3129 


22391808 


22388992 


9.2±2.5 


22±11 


<30 


F2M1532+2415 


22392064 


22389248 


28.4±4.0 


31±13 


<30 


F2M1656+3821 


22392320 


22389504 


9.7±2.4 


15±8 


<30 



Note. — All observations made in cycle 4 (Jul 2007 - May 2008) except for F2M0915-2418, 
whose IRS spectrum was taken in May 2007 as part of program 30121. 
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Figure 1. Stacked IRS rest-frame spectra of red quasars. The 
spectra are normalized to 1 at 15 /im and offset. The dashed lines 
denote the central wavelengths of the family of PAH emission fea- 
tures at 6.2, 7.7, 8.6, 11.2, and 12.7 fim. The shaded region rep- 
resents the typical region (8.0-12.0 fim) where the 9.7 jim silicate 
absorption feature is seen. 

We followed up this sample of 13 red quasars with the 
MIPS and IRS on board the Spitzer Space Telescope (PID 
40143). The IRS observations were made with the Short- 
Low (SL; AA = 5.2-14.5 /im) and Long-Low (LL; AA 
= 14.0-38.0 /im) modules of the IRS, (with some of the 
highest redshift objects not using the SL2 5-7.5 /im mod- 
ule). We also obtained MIPS 24/im, 70/zm and 160/xm 
photometry (with the lowest redshift objects not being 
observed at 160 /Ltm). Table Q] provides the details for 
our Spitzer observations, including the IRS and MIPS 
AORIDs as well as the MIPS 24, 70, and 160 /im fluxes. 

Our infrared spectra and photometry were used for a 



number of purposes. We can infer star formation rates 
from the PAH emission in the IRS spectra and the far- 
infrared excesses from the MIPS photometry. The sil- 
icate feature at 9.7 /im provides us with an estimate 
of the cold dust obscuration, and the mid-infrared con- 
tinuum emission provides an estimate of the intrinsic 
(unreddened) quasar luminosity. Near-in frared and opti- 
cal spectra of these objects are shown in Glikma n et al.l 
(|200ll200l . though we obtained an additional spectrum 
of FTM1507+3129, see Section 4.1.11. 

3. DATA REDUCTION AND ANALYSIS 

The IRS spectra were reduced in the following man- 
ner. First we removed the bad or hot pixels using the 
IDL program irsclean . The resulting spectra were back- 
ground subtracted by differencing the first and second- 
order apertures. They were then extracted and flux cali- 
brated using SPICE version 2.0.1 provided by the Spitzer 
Science Center. Finally, the different orders were stitched 
together using a weighted average in the overlapping or- 
ders, with no significant scaling adjustments needed in 
the order-to-order jumps. A plot of the 13 different nor- 
malized spectra can be seen in Figure [TJ 

We used the post-BCD products for the 24 /im MIPS 
images, as the quality was excellent and there was no 
need to re-mosaic them. For the 70 and 160 /jm data we 
re-mosaiced the data using M0PEX version 18.2.2 with the 
default parameters. We employed two methods for the 
photometric measurements: (a) The APEX package within 
M0PEX with the default parameters and PRFs for point- 
source photometry and (b) SExtractor aperture pho- 
tometry with standard point-source aperture corrections 
provided by the Spitzer Science CenteiQ- We used 10% 
errors for the 70 and 160 /im sources, as M0PEX tended 
to underestimate the errors. 

3.1. IRS spectral analysis 

For all the IRS spectra in our sample we measured the 
equivalent width (EW) of the 6.2 /im PAH emission as 
well as the strength of the 9.7 /im silicate ab sorption fea- 
ture following the methodology presented in lSpoon et al.1 

7 Using the old webpage http://ssc. spitzer. catech.edu/mips/apert 
values can now be fou nd m the MIPS instrument hand book at 
http: / /ssc. spitzer. caltech.edu/mips/mipsinstrumenthandbook/50/ 




Figure 2. The diagnostic diagram of Spoon ct al. (2007) with their sources plotted as small filled circles. Our quasars are plotted as 
larger squares with the PAH values showing arrows being upper limits. Usually sources cluster along two tracks: the diagonal continuum 
track from 1C to 3A (magenta) spanning from Starburst to ULIRG-type systems and the activity track from 1A to 1C (cyan) spanning 
from QSOs to Star burst galaxies . Mos t of the quasars are found in quadrant 1A, but some have Silicate troughs deep enough to be in 
quadrant 2A, which [Spoon et al. (2007) speculates to be transitional objects. 



(|2007l) . In short the PAH EW is obtained by integrating 
above a spline interpolated continuum, while the silicate 
feature strength is inferred from the ratio between the 
observed and continuum flux density: 



Ssii = hi 



f obs (9.7 nm) 
fcont (9.7 [im) 



(1) 



The values are given in Table [5] and plotted in Figure [21 
We discuss the results in the diagram in Section 4. 

3.2. Multiwavelength Modeling 

To obtain a complete view of the physical phenomena 
that are affecting the quasar, we modeled the SEDs of our 
sample from the optical out to the mid-infrared. We used 
the available u, g 7 r : i,z-band SDSS phot ometry^ in the 
optic al. In addition, we used 2MASS fS krutski e et al.l 
l200l) . WISE Band 1+2 (|Wright et al.ll20100 . our MIPS 
photometry and IRS mid-IR spectra. 

8 Except for F2M0729+3336, which has no SDSS coverage. We 
obtained synthetic photometry corresponding to the SDSS filters 
from the optical spectrum 



To model the SE Ds we used a phen omenological model 
based on that of iSaiina et al.l (|2006f ) due to the current 
lack of a good physical understanding of the nature and 
geometry of the dust emitting regions. Many of the pa- 
rameters thus do not have a direct physical interpreta- 
tion. Nevertheless, we considered the effort worthwhile 
as it allows an empirical breakdown of the contributions 
of starlight, hot dust emission from the AGN, and cooler 
dust emission from star formation. The model fits are 
shown in Figure El Full details o f the model compo nents 
as applied to quasars are given in lHiner et al.1 ()2009f ). but 
are summarized below. 

(1) Th e composite quasar spect rum. As the com- 
posite of [Vandcn Berk et all (|2001l) suffers from notice- 
able host galaxy contamination at the long wavelength 
end, we constructed a new compo site by subtracting 
the co ntinuum from the composite of lVan den Berk et al.l 
(|2001[ ). and adding the residual emission line com- 
posite to a new continuum constructed using line-free 
optical/near- in f rare d SED points from the composite of 
I Richards et al.l (|2006f) . In Figure [3] this is the green line. 

(2) A power-law component for the mid-infrared emis- 
sion, with an exponential cutoff at short wavelengths to 
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Figure 3. Broadband rest-frame spectral energy distributions (SEDs) of the red quasars. The five photometry measurements from the 
SDSS are shown as squares, the ones from 2MASS as filled circles, the ones from WISE Band 1+2 as open circles and the ones from MIPS 
as stars (three of the low redshift quasars were not observed at MIPS 160/xm). Also shown is the IRS spectrum with a corresponding grey 
error area. The fitted SEDs of the quasars are shown in red. They are built up from a quasar component (in green) and AGN hot dust 
component (pink) , both of which are reddened by an SMC extinction law. Furthermore a very hot quasar dust component (grey) , a warm 
dust component (cyan) and a cold dust component (blue) are added. The corresponding HST image is shown to the right of each SED as 
an 7 X 7 arcsec grey scale image. 



represent dust sublimation, and a Fermi function cutoff 
at long wavelengths. The functional form used was: 



Lagn = 



,(l-a) e hv/k t 



g^-^hcut)/™ -|- 1 



(2) 



where the normalization, L° AG¥S , power-law index a and 
?sub (a proxy for the sublimation temperature) are al- 
lowed to vary, and t'hcut and w were fixed at 0.11 x 10 14 Hz 
and 0.017 x 10 14 Hz, respectively (h and fcs are the Planck 
and Boltzman constants, respectively). This component 
is represented by the magenta line in Figure [3J Af- 
ter addition, components (1) and (2) were reddened by 
the Galactic Center extinction curve of IChiar fc Tielensl 
(2006) to fit any silicate absorption feature. 

(3) A warm dust component to represent the small 
grain emission from Hn regions, represented by a power- 
law with fixed high and low frequency cutoffs: 



L 



SG 



l sg ^ (1_7) er vlvB%x 



(3) 



this component is poorly constrained in mo st of the fits, 
so 7 w as fixed at a typical value of two (e.g. iSaiina et al.l 
( 2006)). The high and low frequency cutoffs, u sg \ and v sg h 
were set to 0.062 x 10 14 Hz and 0.17 x 10 14 Hz, respectively. 
Figure [3] shows this component in light blue. 

(4) a thermal greybody model for the far-infrared emis- 
sion: 



J FIR 



FIR ' 



,3+0 



e hv/k B T _ l 



(4) 



For all of the objects, we fixed the temperature of this 
component, T, to 40K, because of the uncertainty in 
our 160/zm photomet ry. /3 was fixed at 2.0 as above 
(jDunne fc Ealcs 2001). In Figure |3J this is represented 
with a dark blue line. 

(5) Many quasars objects also required an extra near- 
infrared component to match the SEDs. This "very hot" 
component was modeled as a 1000K black body (equa- 
tion 4 with T = 1000). The physical origin of this com- 
ponent is unclear, but it is most likely due to hot dust 
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Figure 3. Cont. 



close to the sublimation radius (e.g..lGlikman et al.ll2006t 
INetzer et all [20071: IMor fc Trakhtenbrot 1 1201 ID . This 
last component can be seen in grey color in Figured] 

In addition we trie d to add a PAH model as described 
in lLacv et all (|2007bl ) , however PAH emission was rela- 
tively minimal and did not improve our fits. We there- 
fore excluded this component in the fitting process. The 
contribution from starlight is so small compared to the 
quasar and star formation components at all wavelengths 
except bluewards of «5000A, so we decided to ignore the 
stellar contribution in our fitting process. 

Most of our quasars are well-fit by this model, with 
relatively small X /dof values, but in some cases the 
fit at the near-IR regions of the spectrum described by 
2MASS magnitudes provides is poor. This is most likely 
due to the host galaxy emission dominating the optical 
and/or near- infrared passbands. 

In the cases where the model fit the broadband SED 
well, the reddening values derived from the Silicate ab- 
sorption feature were in general agreement to the val- 
ues of E(B — V) deri ved from the SMC continuum fit to 
the o ptical spectrum ([Glikman et al.|[2007l lUrrutia et al.l 
2009) , confirming cold dust as the absorbing and redden- 
ing mechanism in these systems (see Table [2]). How- 
ever, in four cases (F2M0915+2418, F2M1151+5359, 
F2M1507+3129 and F2M1532+2415) there was no Sil- 
icate in absorption to account for the reddening seen in 
the optical and near-infrarcd. In those cases, we artifi- 
cially fixed the extinction ry that would've been derived 
from the fit from the IRS spectrum manually. These ob- 
jects have their modeled extinctions marked with a star 
in Table [5] and their fixed extinctions labeld in Figure [3J 

We then derived the obscured quasar luminosities by 
integrating over the quasar components (l) + (2); the in- 
trinsic luminosities quoted in Table [2] are then derived by 



accounting for the extinction obtained from the Silicate 
absorption fit. Since the AGN contribution is well con- 
strained by the IRS spectrum and optical/near- IR pho- 
tometry points, the uncertainties are overall small (< 0.2 
dex) and only pertain to the accuracy of the Silicate ab- 
sorption feature. The FIR/Starburst contribution was 
taken solely from the integration of the two dust bumps 
(light and dark blue; (3)+(4)), not from any component 
in the optical regime. Since that FIR contribution is 
largely determined by the 70 and 160 /im photometry 
points alone, it has the largest uncertainties in the contri- 
bution (~ 0.5 dex). The Starburst contribution ranges 
from log(LFiR,/£o) = 10.67 - 12.67, placing the host 
galaxies of red quasars roughly in the LIRG regime. For 
the QSO contribution, we find them to have bolomet- 
ric luminosities in the range of log(LQso/-£<©) = 12.03 - 
13.37. These are very large values, placing these quasars 
among the brightest in the Universe in the IR further- 
more strengthening our previous conclusion that we are 
only probing the tip of the quasar luminosity iceberg 
(jGlikman et al.l[2004 [20071 lUrrutia et al.ll2009t i. 

3.3. Black hole masses and accretion rates 

We estimated black hole masses (Mbh) for the twelve 
of our thirteen quasars for which broad lines were visible 
in either the optical or near-infrared spectra (Table [3]) . 
H/3 was our preferred line as it is relatively isolated, and 
most quasar black hole mass estimates are calibrated on 
it, but in most cases the broad H/3 line was too heavily 
extinguished, so a longer wavelength line such as Ha or 
Pa/3 was used (Pa/3 was generally preferred as it is free 
of contaminating narrow lines). This assumed that the 
H/3, Ha and Pa/3 lines have the same velocity widths, but 
was indeed the case (within the errors) for the objects in 
which we could measure two or more broad line widths 
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Table 2 

Parameter results from multiwavelength fitting 



Source 


QSO contrib- 


Starburst (FIR) 


QSO / SB 


E(B-V) a 


E(B-V) b 




D c 




ution (intrinsic) 


contribution 


Ratio 


(spectrum) 


(SED fit) 




coefficient 




(log Lq SO /Lq) 


(log L fir /Lq) 












F2M0729+3336 


12.83 


<11.47 


> 22.91 


0.83±0.22 


0.605 


-0.43± 0.09 


1.92 


F2M0825+4716 


13.10 


12.35 


5.62 


0.52±0.10 


0.806 


-0.50± 0.03 


1.29 


F2M0830+3759 


12.03 


11.35 


4.79 


0.80±0.15 


0.478 


0.00± 0.06 


1.39 


F2M0834+3506 


12.04 


11.42 


4.17 


0.58±0.05 


0.629 


0.14±0.03 


0.95 


F2M0841+3604 


12.12 


12.03 


1.23 


1.34±0.11 


0.770 


-0.76± 0.15 


1.73 


F2M0915-2418 


13.37 


12.69 


4.79 


0.36±0.12 


0.354* 


-0.05± 0.05 


1.29 


F2M1012+2825 


12.24 


<11.88 


>2.29 


0.82±0.10 


0.647 


0.03± 0.10 


1.13 


F2M1113+1244 


13.14 


12.43 


5.13 


1.41±0.11 


0.294 


-0.31± 0.01 


1.24 


F2M1118-0033 


12.71 


12.29 


2.63 


0.85±0.11 


0.512 


-0.15± 0.04 


1.40 


F2M1151+5359 


12.15 


11.77 


2.40 


0.42±0.08 


0.532* 


-0.01± 0.07 


0.95 


F2M1507+3129 


12.57 


12.05 


3.31 


0.74±0.09 


0.532* 


-0.01± 0.07 


1.16 


F2M1532+2415 


12.47 


10.67 


63.10 


0.90±0.53 


0.354* 


-0.09± 0.02 


1.50 


F2M1656+3821 


12.44 


10.83 


40.74 


0.88±0.16 


1.169 


-1.31± 0.01 


1.34 



a Derived via spectral fitting of a quasar composite reddened by the SMC law, either from Urrutia ct al. (2008) or Urrutia ct al. 
(2009); the latter involves much more careful fitting, masking emission line regions. 

^ Derived from the SED fitting and modeling as described in Section l3.2l the extinction (tv) is derived from Silicate absorption, 
except for objects marked with a star. Note the general agreement with the Silicate depth derived via the method of Spoon ct al. 
f200% . 

c D — G/C, used as disturbance parameter. Values of G and C arc from lUrrutia et alj 12003) , sec that publication for details 
on how the values were obtained. D values close to 1.0 represent morphologically undisturbed systems. The higher the D value, 
the more disturbed the system is. 



(F2M0825+4716, Ha and Pa^; F2M0834+3506, H/3 and 
Ha; F2M0915+2418, Ha and Pa/3; F2M1113+1244, H/3, 
Ha and Pa/3). 



Black hole mass estimates were then calculated using 
the broad-line widths and quasar luminosity assuming a 
broad line region radius, Rblr- 

Rblk = K + alog(AL A (5100 J 4)), (5) 

(|Bentz et all 120091) with K = -21.3 and a = 
0.519 and estimating the intrinsic luminosity at 
5500A ,AL A (5100,4) by measuring the quasar luminos- 
ity at a rest-frame wavelength of 15^m (where extinction 
should be negligible) and applying the ratio of the bolo- 
metric corrections at 5100A (12) to the correct ion at 
rest-frame 15^m (9.0) from (jRichards et al.ll2006t ). The 
black hole mass estimate, A/bh is then given by: 

M BH = f 7 £~ (6) 

where a is the velocity dispersion of the broad line, G is 
the gravitational constant, and we assume a geometric 
correction factor of / = 5.5. The Eddington luminosity 
is: 

L Edd = 1.3xl0 31 (^)w, (7) 

and Eddington ratios are estimated by dividing the in- 
trinsic luminosity of the quasar from the fit by the Ed- 
dington luminosity. 

4. RESULTS 

We first place our red quasars into th e context of 
other luminous infrared galaxies using the iSpoon et al.1 
(2007) diagnostic diagram (Figure [5]). This plot of sili- 
cate absorption depth versus PAH equivalent width sep- 
arates the known classes of infrared-luminous galaxies 



into two. One sequence extends from region 1C (com- 
prised mostly of starburst dominated objects) into 1A 
(comprised mostly of IR-bright quasars) as continuum 
replaces PAH emission in objects with little silicate emis- 
sion or absorption. We plot this sequence with cyan dou- 
ble arrows in Figure [5] The second sequence extends 
from 1C through 2B through 3A in which objects again 
become continuum dominated, but with high silicate ab- 
sorption depths (the "warm ULIRG" population) , shown 
with magenta double arrows in Figure [21 

Normal quasars are typically found in region 1A, as 
are the majority of our dust-reddened objects, but the 
fact that we see one object firmly and one barely in the 
the sparsely populated quadrant 2A of the diagram is 
noteworthy. Such objects are rare in most compilations 
of galaxies and AGN and are believed to be part of a 
transitioning population, confirming that we have indeed 
located an interesting class of object, perhaps intermedi- 
ate between the heavily absorbed warm ULIRGs and the 
normal quasar population. Also, it is worth mentioning 
that none of the sources show silicate in emission, which 
is unusual for these very IR-bright quasars; it's possible 
that in some of the objects (like the four we had to fix 
the Silicate absorption manually in the SED fitting) the 
silicate absortion is missing because it is being filled with 
an emission bump. 

The results of the fitting (Figure [3] and Table [2]) show 
that the AGN component dominates the bolometric lu- 
minosity in all cases, with most of our objects having 
measurements or limits on their far-infrared luminosi- 
ties corresponding to star formation rates of ~ 10 — 
SOOMqy^ 1 using the iKennicutt et al.l (|1998| ) conversion 
between far-infrared luminosity and star formation rate. 
Such star formation rates indicate luminous starbursts, 
but not at the level at which they are adding significantly 
to the stellar mass of the host galaxies. 

Next we consider the accretion luminosities of our ob- 
jects. Our objects span a wide range of luminosities rela- 
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Table 3 

Black hole mass estimates and Eddington rates 



Object 


Redshift 


Line 


FWHM 


log 


log 


log 


logio 


logio 










(km s —1 ) 


(erg s _1 ) a 


(Mbh/M q ) 


(L/ I/Bdd) 


(Lb/Lq) 
(fit) b 


(Lb/Lq) 
(monotonic) c 


F2M0729- 


s-3336 


0.957 


H/3 


28oo± 200 


45.8 


8.7±0.1 


O.liO.l 


11.06 


10.99 


F2M0825- 


^4716 


U.oUO 




ZOO*±nz 1UU 


AP, 1 
^40. 1 


O.OIEU. ± 




1U.OU 


ru.yu 


F2M0830- 


[-3759 


0.414 


H/3 


3741± 200 


45.3 


8.6±0.1 


-0.4±0.1 


10.97 


11.01 


F2M0834- 


s-3506 


0.470 


Ha 


10500±2000 


45.2 


9.5±0.1 


-1.3±0.1 


10.69 


10.65 


F2M0841- 


^3604 


0.553 


Ha 


3380± 500 


45.2 


8.5±0.1 


-0.4±0.1 


10.53 


10.73 


F2M0915- 


^2418 


0.842 


Pa/3 


4186±800 


46.5 


9.4±0.2 


0.1±0.1 


10.58 


10.66 


F2M1012- 


^2825 


0.937 


Ha 


8050± 1000 


45.3 


9.3±0.1 


-1.0±0.1 


10.62 


10.79 


F2M1113- 


1-1244 


0.681 


H/3 


2276± 300 


46.1 


8.6±0.1 


0.4±0.1 


10.91 


10.91 


F2M1118- 


-0033 


0.686 


Pa/3 




45.7 


<8.7 d 




11.01 


11.12 


F2M1151- 


1-5359 


0.780 


Ha,Pa/3 


5700±500 


45.4 


9.0±0.1 


-0.7±0.1 


10.48 


10.88 


F2M1507- 


1-3129 


0.988 


Ha 


12800±2000 


45.7 


9.9±0.1 


-1.3±0.1 


10.60 


10.96 


F2M1532- 


^2415 


0.562 


Pa/3 


11700±1000 


45.6 


9.8±0.1 


-1.2±0.1 


10.95 


11.09 


F2M1656- 


^3821 


0.732 


Pa/3 


3064±300 


45.4 


8.5±0.1 


-0.2±0.1 


10.54 


10.67 



a monochromatic luminosity at 15/xm (rest). 

rcst-framc S-band luminosity of the central component of the host galaxy as given in table 6 of lUrrutia et al.l 1 120081) ("Total" 
host magnitude in table 6). 

c rest-frame S-band luminosity of the extended host galaxy (all components) ("mono" luminosity in table 6). The B-band 
luminosities were estimated from the observed /-band fluxes assuming a stellar SED flat in F\ to perform the small k-corrcctions 
required to map observed /-band to rest-frame £?-band. 

limit on the black hole mass assuming Eddington-limited accretion. 



tive to the Eddington Limit (^Edd), including many ob- 
jects with high luminosities when c ompared to typical 
quasa r Eddington Ratios ~ 0.1 (e.g., Stcinhar dt fc Elvisl 
2010). This could be an artifact of partial obscura- 
tion, for example, if the inner broad line region, where 
the linewidths are greatest, remains somewhat obscured. 
However, where we are able to measure the widths of 
multiple broad lines in the same object, the linewidths 
are very similar irrespective of wavelength, making this 
latter possibility unlikely. 

Figure Q] displays the black hole mass of our quasars 
versus their host luminosity, left for the central bulge 
component and right for the total host light. The cen- 
tral bulge luminosity was derived by fitting the bright- 
est host galaxy components, while the total luminosity 
component included the flux in disturbed, outlying com- 
pon ents, both aft e r sub traction of the quasar nucleus 
(see lUrrutia et al.l (|2008T ) for an in depth discussion on 
the AGN - host galaxy fitting). Objects with high ac- 
cretion rates (above 0.3 Eddington) are plotted in red, 
objects with low rates in blue. The Figure shows that 
there is very little correlation between the black hole 
masses and the corresponding galaxy luminosities, either 
of the fitted central galaxy component or to the total host 
galaxy light. However, none of our black hole masses 
exceed the black hole mass - bulge luminosity relation 
(M arconi fc Hunt 2003) by a larg e factor, also when cor - 
rected for evolution according to IBennert et al.l (|2010f ). 
even when only the central component of a merging sys- 
tem (defined as the one containing the quasar nucleus) 
is considered. The fact that we find a significant fraction 
of our objects with high accretion rates below the black 
hole mass - bulge luminosity relation leads us to specu- 
late that they are young, just ignited black holes. Their 
growth has started after the host galaxy growth via a 
starburst event, yet since they are accreting at such high 
rates they will move onto the "classical" relation within 
~10 7 years. 

In Figure [5] we plot both the total FIR-luminosity 



against the AGN luminosity, and the ratio of accretion 
rate to the Eddington limit of the quasar against the 
ratio of the far-infrared luminosity to the £?-band lumi- 
nosity of the host galaxy (a proxy for the specific star 
formation rate). There is a hint of a correlation between 
the total FIR lumi nosity and AGN lum inosity (as found 
for nearby quasars. iNetzer et al.ll200"7l ). though it is not 
tight. There seems to be no relation between the mass- 
normalized quantities of specific star formation rate and 
accretion rate relative to the Eddington Limit. 

Finally, in Figure |6l we examine the relationship be- 
tween the host galaxy morphology and the star formation 
rate and AGN accretion rate. In order to quantify the 
degree of disturbance of the host we use the ratio of Gini 
C oefficient (G) t o Conc entration Index (C). We showed 
in lUrrutia et al.l (|2008f) that the combination of these 
two quantities gave a robust indication of the degree of 
disturbance of the host galaxy. We therefore define a 
"disturbance parameter", D = G/C. This quantity is 
appro ximately unity for normal galaxies ( Abraha m et all 
2003), but is significantly greater than unity for highly 
disturbed systems. There is a weak (though not statis- 
tically significant) trend for high accretion rate objects 
to have more disturbed morphologies, though there is no 
apparent trend with specific star formation rate. 

5. DISCUSSION 

Our aim in this study was to determine whether the 
infrared properties of these dust-reddened objects were 
consistent with the merger-driven co-evolutionary model 
for the growth of SMBH and their host galaxies. In 
particular, whether these objects are an intermediate 
stage between a merger-induced starburst and a "nor- 
mal" quasar with an unreddened, blue continuum, with 
AGN feedback having stopped further star formation in 
the host. Our results seem to be consistent with this 
overall picture, but we see considerable variation in prop- 
erties such as silicate absorption depth, star formation 
rate and AGN accretion rate on an object-by-object ba- 



Spitzer observations of young red quasars 



9 





10,4 10.6 10.8 11.0 11.2 
\og m {L B /L Q ) (central host component) 




10.4 10.6 10.8 11.0 11.2 
log 10 (L B /L G ) (total host luminosity) 



Figure 4. Black hole mass versus host luminosity for our red quasars. Left: using the model fits to the brightest host galaxy component of 
lUrrutia et al. (2008), right: using the total host luminosity, including the flux in disturbed, outlying components not well-fit by the model. 
Cyan symbols denote objects accreting at <30% of the Eddington rate, red points objects with accretion rates >30% Eddington, the object 
with a limit only on its black hole mass (assuming Eddington-li mited accretion ), F2M1118-0033, is shown as a black dot. The red dashed 
line is the local bl ack hole mass — bulg e luminosity relation from Ivan dcr Marcl (1999), the black dashed line is the relationship at z ~ 0.7 
(interpolated from Merloni ct al. 2010). Most of our high accreting object lie below the black hole mass - bulge luminosity relation. 



13.0 



12.5 



12.0 



2.5r 



11.5- 



11.0 
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Figure 5. The log of the ratio of far-infrared luminosity to host luminosity in B-band (as a proxy for the specific star formation rate) is 
plotted against the accretion rate relative to the Eddington limit. Symbols are as in Figure 4. While the most luminous quasars also show 
the highest luminosity in the FIR, there is no correlation between the inferred star formation and accretion rate of the quasar. 

"young quasars" . Not all dust-obscured quasars fit in the 
co-evolutionary scenario, but statistically, it is the most 
likely explanation and therefore our conclusions are to 
be taken as such. 

All these objects have mid-infrared (AGN) luminosi- 
ties in excess of their far-infrared (presumed starburst) 
luminosities. This dominance of the AGN suggests that 
AGN feedback processes, if they are indeed important, 
will have started to have their impact on star forma- 
tion in the host galaxy as we would have expected 
the starburst component to have dominated otherwise. 
We do not see any obvious quasar feedback "smok- 
ing guns"; broad absorption line features from AGN 



sis. 

Much of the object to object variation is due to the 
heterogeneity of our sample; some have high reddenings, 
morphological disturbances, some do not. While they 
all are luminous quasars, they do span a large range in 
luminosities, too. In the overall picture, however, trends 
are evident. If an object is morphologically disturbed, it 
is much more likely to have a high Eddington ratio and to 
show Silicate in absorption (be reddened by cold dust). 
Similarly, the two quas ars which did not sho w any signs 
of merger interaction in lUrrutia et al.l (|2008[) also do not 
display high Eddington ratios or Silicate in absorption 
and therefore should not be included characterized as 
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Figure 6. Right: Host galaxy morphological distortion, parameterized as D, the ratio of the Gini Coefficient (G) to the Concentration 
Index (C) against Eddington ratio. Normal galaxies have a mean D 1.2 (Abraham et al. 2003), objects with values below this (shaded in 
the plot) are relatively undisturbed, those with higher values have increasing amounts of disturbance. Symbol colors as for Figure [4] Left: 
Morphological distortion versus the log of the ratio of far-infrared luminosity to host luminosity in B-band, an indicator for star formation. 
While there could be a hint of a correlation between the amount of host disturbance and the accretion rate, there is no correlation found 
for star formation versus the host disturbance. 

sample, J1356+1026 , is a relatively strong radio source 
(jGreene et al.l l2011bf ). Clearly the role of radio jets in 
powering these outflows needs to be established. 

In agreement with the evolutionary model, we see a 
large fraction of objects with high black hole accretion 
rates, consistent with the idea that quasar accretion rates 
are high in the early phase of growth. The objects with 
the highest accretion rates also tend to fall below the 
local black hole mass - bulge luminosity relation by ~ 
0.5dex; none of our objects with lower accretion rates 
are found there. This offset from the local relation is 
consistent with a picture in which the black hole is yet 
to grow to its equilibrium size following a major merger. 
It is also broadly consistent with numerical simulations 
of black hole growth in quasar merger scenarios, in which 
the total black mass hole gro ws by about 0. 5 dex with 
every major merge r event (e.g.. lLi et al.ll2007l figure 10). 
iSarria et al.l (|2010l ) also find that their luminous, dust- 
obscured quasars at z ~ 2 fall below the black hole mass 
- bulge mass relation at that redshift. An interesting 
exception is F2M0915+2418, which has a high accretion 
rate, but falls close to the black hole mass - bulge mass 
relation. This object also has one of the highest star 
formation rates in the sample, however, allowing it to 
evolve along a trajectory that will keep it close to the 
black hole mass - bulge mass relation. 

We do, however, have some remaining puzzles, where 
our data do not fit in with t he most naive e xpect ations of 
the evolutionary model. In lUrrutia et al.l ()2008| ) we saw 
a relationship between D and the amount of reddening in 
the host, but in Figure [6] there seems to be little correla- 
tion between the degree of disturbance of the host galaxy 
and the star formation rate, except that all high accreting 
quasars are also disturbed to some degree. There may, 
however, be a very weak trend for objects with high ac- 
cretion rates relative to the Eddington Limit to have high 
values of D, see Figure|6l Similarly, there is no clear mor- 



thermal winds, for example, are not present in any of 
our objects which are of high enough redshift to see 
Mg II absorp t ion ( z ^0.7). However, experience from 
lUrrutia et al.l (|2009l ) , has shown that for the red quasars 
selected using our criteria the spectrum below w 5000 
A is so extinguished that even deep absorption features 
are difficult to detect in objects with z <0.9. Only 
three of our objects (F2M0729+3336, F2M1012+2825 
and F2M1507+3129) have appropriate redshifts, and 
F2M0729+3336 is too noisy to detect absorption fea- 
tures. Nevertheless, an excess of low-ionization, broad 
absorption line qu asars has been found in our full sam- 
ple of red quasars. lUrrutia et al.l (|2009l ). and a tentative 
relationship between wind properties and star formation 
seen has been seen bv lFarrah et al.l (|2012l ) in samples of 
reddened quasars. 

Blueshifted wings to [O III]5007A, with veloc- 
ities relative to the line peak of up to ~1000 
km s^ 1 are seen in five of our thirteen ob- 
jects (F2M0729+3336, F2M0825+4716, F2M0915+2418, 
F2M1118-0033, F2M1151+5359). Of the four of these 
for which we can estimate accretion rates, three have ac- 
cretion rates just exceeding the Eddington limit, possi- 
bly indicating a link with extremely high accretion rates, 
but the star formation rates in the hosts vary widely. We 
also lack any information on the spatial extent of these 
outflows, which may be confined to the nuclear regions. 
Blueshifted [O III] line components a re common in com- 
pact radio sources: iHolt et aLl (|2008t ) propose that these 
components are due to outflows of gas driven by the radio 
jets. The radio sources in our quasars are however sig- 
nificantly l ess luminous than those in the radio galaxies 
studied bv Holt et all ([2008ft . Radio-quiet type-2 lumi- 
nous quasars at low reds hift also tend to show outflows 
m [O III]5007A emission (|Greene et al.ll2011al) . suggest- 
ing that these are not confined to radio sources, though 
the most spectacular object in the lGreene et al.l (|2011a[ ) 
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phological trend related to the offset of the quasar from 
the local galaxy luminosity - black hole mass relation. 

To some extent, the lack of clear trends here may 
be due to an admixture of objects reddened by fore- 
ground galaxies. In particular the low accretion rate 
objects may include some objects where the reddening 
is from an intervening galaxy through a chance align- 
ment, perhaps from a galaxy in the same group or cluster. 
F2M1532+2415 is a good example of this, with an appar- 
ently highly-disturbed, early stage merger morphology, 
but both a low accretion rate and a low star formation 
rate, and F2M0834+3506 is a very good candidate for a 
low accretion rate broad-line radio galaxy reddened by a 
foreground irregular galaxy. 

In summary, our observations are consistent with the 
simplest evolutionary models, where a merger induced 
starburst triggers a obscured quasar which begins life ac- 
creting at close to the Eddington rate and later evolves 
into an unobscured quasar accreting at the ~ 0.1 of the 
Eddington rate more typically observed in the normal 
quasar population. We see no obvious evidence, how- 
ever, of feedback affecting star formation in this particu- 
lar sample (though we do see evidence of strong outflows 
of ionized gas in some objects). Nor do we see a clear se- 
quence or progression from a newborn, highly obscured, 
high accretion rate AGN with a very disturbed host into 
a quiescent, almost unreddened object with a lower ac- 
cretion rate. Our quasars were, however, selected within 
a fairly narrow range of luminosity (due to the 2MASS 
detection requirement and restricted redshift range) and 
reddening (if the reddening were too high, they would be 
classed as type-2s, but if it were too low they would not 
be picked out as red quasars), so we might not expect to 
be able to see much of an evolutionary sequence. 

These observations are most consistent with a pic- 
ture in which the majority of objects undergo the star- 
burst/merger phase well before we see these reddened 
quasars emerge, though not so long after that evidence 
for morphological disturbance of the hosts has been 
erased. The evolution in the black hole mass - bulge 
mass plane is then a shift along the bulge mass axis by 
~ +0.5dex followed ~ 10 8 yr later by a shift along the 
black hole mass axis by a similar amount. In one case, 
however, (F2M0915+2418) we do seem to see the star- 
burst and accretion happening simultaneously. As our 
selection techniques for the dust-reddened quasar popu- 
lation expand, particularly as mid-infrared selected sam- 
ples become large enough to include significant numbers 
of reddened type-1 quasars we expect to be able to ex- 
pand our samples considerably, and trace out more of the 
evolutionary scenario. 
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APPENDIX 

APPENDIX: NOTES ON INDIVIDUAL OBJECTS 
F2M0729+3336 

This object lies outside the SDSS coverage, so photometry from the ESI spectrum was used when constructing the 
SED. Both the silicate feature and the optical SED are well fit with a reddening of E(B — V) = 0.83 ±0.22. The object 
is detected in the X-ray by Chandra with six hard X-ray photons (jUrrutia et al .1120051 ). It has an accretion rate close 
to (formally slightly in excess of) the Eddington limit, but no evidence of cool dust from substantial star formation 
activity, though the host galaxy appears disturbed. The inferred black hole mass lies significantly below the black hole 
mass - bulge luminosity relation, consistent with our observation that the black hole is still growing rapidly. 

F2M0825+4716 

Another high accretion rate object, but this one also has a strong starburst component in the far-infrared. The host 
is very disturbed, but, if the fit to the central component of the galaxy only is considered, lies close to the black hole 
mass - bulge luminosity relation. The strong black hole growth would thus take it significantly above the relation 
within a Salpeter time, if it were not for the high star formation rate, and the likelihood that some of the extended 
components of the host galaxy would merge with the central one in that timescale. As such, F2M0825+4816 is a prime 
example of black hole - host galaxy co-evolution. 

F2M0830+3759 

The lowest redshift (z=0.414) object in our sample, with a strong Chandra detection (|Urrutia et a l. 2005). A deeper 
X-ray spectrum with XMM reveals a cold or only moderately ionized absorber with a column density consistent with 
it being ass ociated with the dust causing the quasar reddening, as well as an Fe Ka emission line from a reflection 
component (jPiconcelli et al.ll2010D . This object has a relatively massive host galaxy, and the object lies below the 
black hole mass - bulge luminosity relation. Nevertheless, the accretion rate is relatively high L/L^dd ~ 0.4, so again, 
this quasar may be able to build black hole mass effectively. The star formation rate is constrained to be relatively 
low, ~ lOMoyr" 1 . 

F2M0834+3506 

This object has little or no silicate absorption, inconsistent with the observed reddening in the optical. The HST/ACS 
image too suggests that this, rather than being a post-merger system, may be a normal quasar obscured by a foreground 
irregular galaxy. The spectrum shows multiple-peaked permitted lines, with three identifiable components, one at the 
same redshift as the forbidden lines, one at w +3600kms _1 and one at « — 5100kms _1 . The object is the only one 
of two i n our sample to show mult iple-peaked permitted lines, typically only seen in low accretion rate, radio-loud 
objects (jEracleous fc Halpernl 120031 ) . and our analysis indeed confirms a low accretion rate and a relatively low star 
formation rate. It lies a little above the local black hole mass - bulge luminosity relation, consistent with the idea that 
this object is a fairly normal quasar obscured by dust from a foreground object. 

F2M084 1+3604 

Consistent the very irregular morphology, this object is well-fit by the simple cold dust screen model. It has both 
a moderately-high accretion rate and moderately-high star formation rate, lying below the black hole mass - bulge 
luminosity relation. The radio informati on reveal a steep spectrum point source, while the X-ray image shows a very 
hard hardness ratio (jUrrutia et al.ll2005h . 

F2M0915+2418 

This object is the most luminous quasar in our sample, and is bright enough in the i nfrare d to be detected in the 
IRAS Faint Source Catalog. The IRS spectrum is an archival one from ISargsvan et~a l. (2008). This object has both 
a high star formation rate and high accretion rate, and is the only one of our high accretion rate objects to fall close 
to the black hole mass - bulge mass relation in both plots of Figure |4] Due to its high star formation rate and high 
accretion rate, it will evolve with a trajectory close to the black hole mass - bulge mass relation for at least a Salpeter 
time (~ 10 8 yr). 

F2M1012+2825 

A fairly unexceptional object in terms of its accretion rate, star formation rate and position on the bl ack hole mass 
- bul ge luminosity relation, it does, however have a double nucleus with a separation of only 1.2kpc (jUrrutia et al.l 
f2008h . 

F2M1113+1244 

Another good candidate for a young quasar, with a high accretion rate (super Eddington; log(L/ L^dd) = 0.4 ±0.1), 
high star formation rate (FIR contribution in ULIRG regime; log(Lpm/ Lq) — 12.39), and plotting well below the 
black hole mass - bulge mass relation. 
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Figure 7. Near-infrared spectrum of F2M1507+3129, showing a double-peaked H« line. 

F2M1 118-0033 

This object shows no clear broad lines in either the optical or near-infrared spectrum. Ha is detected, but is relatively 
narrow (« lOOOkms -1 ). There is a hint of broad Pa/? in the near-infrared if-band, but it is not strong enough to allow 
us to estimate a linewidth. The SED is consistent with a reddened quasar. No black hole mass estimate is possible in 
this case, though we have obtained a limit based on the quasar luminosity and assuming Eddington-limited accretion. 

F2M1151+5359 

The reddening of this quasar is relatively modest, and the host almost undisturbed. Nevertheless, it has a high 
accretion rate and relatively high star formation rate, and lies well below the black hole mass - bulge mass relation. 

F2M1507+3129 

This object has a near- infrared spectru m inlGlikman et all (|2004j ). however, we obtained an improved near- infrared 
spectrum with the Triplespec instrument (jHerter et al.H20d8f) on the Hale 200- inch telescope on July 11th 2011. The 
data were calibrated using the AOV s tar HD127304. The new spectrum clearly shows a double-peaked broad Ha 
emission line in J-band (Figure lA.llf) . The object is the second highest redshift double-peaked quasar currently 
known, and is m ore luminous than CXOECDFS J0331-2755 at z = 1-369, the most distant double-peaked emitter 
(|Luo et al.ll2009l) . Like many objects in this class (|Lewis fc Eracleous Il200d ). and indeed similar to the other double- 
peaked emitter in our sample, F2M0834+3506, it has a relatively low accretion rate (0.05 Eddington). 

F2M1532+2415 

The SED of this object is difficult to fit, consistent with the complicated structure seen in the host galaxy image. 
Despite a red continuum, the IRS spectrum lacks silicate absorption. The far-infrared luminosity is very low, the 
accretion rate is also low and the object lies close to the black hole mass - bulge mass relation. As of now the 
underlying physics responsible for the reddening and host disturbance is difficult to explain. 

F2M 1656 +3821 

This object has a relatively low star formation rate (barely above the LIRG regime; Iog(LFiR/£©) = 11-19), but a 
high accretion rate (roughly Eddington rate), and should be able to accrete matter onto the black hole to move up 
to the black hole mass - bul ge mass relat i on fro m its current position w 0.5 dex below. It is also the object that is 
clearly in quadrant 2A of the iSpoon et al.l (|2007| ) diagram, since it has such a deep Silicate absorption feature. 



